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^ Abstract. Wc study transport length scales in carbon nanotubes and graphene ribbons under 

fH the influence of Anderson disorder. We present generalized analytical expressions for the density 

I of states, the elastic mean free path and the localization length in arbitrarily structured quantum 

<yi wires. These allow us to analyze the electrical response over the full energy range, including the 

mJ regions around van Hove singularies, traditionally difficult to access by alternative approaches. 

M Gomparing with the results of numerical simulations, we demonstrate that both the diffusive and 

the localized regime are well represented by the analytical approximations over a wide range of the 
energy spectrum. The approach works well for both metallic and semiconducting nanotubes and 
nanoribbons but breaks down near the edge states of zigzag ribbons. 
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The exceptionally high electrical conductivity of carbon nanotubes (CNTs) can be attributed to a 
combination of various factors. The quasi-one-dimensional (quasi-lD) crystalline structure allows the 
preparation of samples with a very low defect density [1]. The rigidity of the sp^-hybridized carbon 
lattice minimizes the effects of thermal vibrations [2]. One truely unique feature of metallic CNTs, 
however, is their electronic structure [3]: independent of the tube diameter, there are exactly two 
massless bands of high velocity crossing at the Fermi energy, resulting in a very low density of states 
(DOS) which effectively suppresses scattering even in the presence of disorder [4] and allows ballistic 
transport over hundreds of nanometers [5, 6]. 

Planar, finite-width graphene nanoribbons (GNRs) display many similarties to their roUed-up 
counterparts. The band gap of armchair-edge GNRs depends on the width just as it depends on the 
circumference in zigzag CNTs. Zigzag-edge GNRs, on the other hand, are always metallic just as 
armchair CNTs. A remarkable difference, however, lies in the edge states generally found at zigzag 
edges in graphene [7, 8] that are nearly localized by their extremely low dispersion and exist beside the 
conducting channels at the band center [9]. Energetically, these states lie exactly around the Fermi 
energy, so it would be expected to have an important impact on the low-energy transport properties 
of GNRs. The low velocity of the edge channels results in a large DOS, which should by itself lead to 
strong scattering. However, the local DOS (LDOS) is strongly concentrated on the edge atoms and 
protecting the conduction channel from scattering into the edge states. 

The theory of quantum transport in disordered systems has been studied for many decades [10]. 
Typically, the underlying model for such studies is that of free electrons or, if a discretization is 
desired, electrons on a (square) lattice. Quantities such as effective electron mass, density of states or 
Fermi velocity are then viewed as free parameters that can be adjusted to fit the properties of metals 
[11]. Graphene and CNTs, however, have an electronic structure that is very different from any 
effective-mass free-electron approximation. In two-dimensional graphene, the electrons have a Dirac- 
like massless dispersion [12, 13], which results in equally massless bands of high velocity in metallic 
CNTs [3] . Tuning the Fermi energy away from the charge neutrality point (CNP) via electronic gating 
or chemical doping, CNTs of all chiralities go through a sequence of additional bands that lead to 
van-Hove singularities in the DOS which are typically not covered by simple models based on a smooth 
or fixed DOS and number of conductance channels. 

In this paper, we will present a generalized expression for the DOS in arbitrarily structured, 
Anderson-disordered quantum wires that allows to access the whole energy range including the regions 
around the van Hove singularities. Using this representation of the DOS, we then show that a general 
expression, originally devised for metallic quantum wires with a fixed number of channels still holds 
true near band edges if used correctly. The various expressions will be applied to both CNTs and 
GNRs, including the special case of pure edge disorder in the latter case. 

For ID systems such as CNTs, it was argued by Mott and Twose [11] that all quantum states 
should become localized at arbitrarily low disorder, leading to an exponential suppression of the 
conductance at zero-temperature. For truly ID quantum wires, Thouless demonstrated that the 
localization length £ioc is identical to the mean free path H^i [14]. Later, he derived a generalized 
expression for metallic quantum wires of finite diameter [15], 

Aoc ~ ^4l, 1 

37r^ 
where the prefactor consisting of cross section A and Fermi wave vector fcp can be identified as the 
number of conduction channels A'ch in the modern language of quantum transport theory. A more 
accurate expression derived by Beenakker from random matrix theory reads [16] 

Aoc« [/3(A^ch-l)/2+l]4i, (2) 

where /3 = 1 for time-reversal invariant systems and (3 — 2 otherwise. In either case, however. Id and 
A^ch are viewed as fixed parameters and it is not at all clear, how these relations should be applied to 
systems where these quantities vary strongly with the energy and may not even be well defined near 
band edges where N^h is discontinuous. An alternative approach based on a perturbative expression 
for the Lyapunov exponents takes into account the true band structure of a quantum wire, but it still 
fails to reproduce the correct behavior near band edges [17]. 



Diffusion and localization in carbon nanotubes and graphene nanoribbons 



& 

^ 



(5,5) CNT 

a^ =0.5 eV 


II ! 1 
II 1 1 




Z^"*^-*--— — _— /' ^^"^ 




J^ i 


1 






E{eY) 



Figure 1. Density of states g{E) and localization length li^^ of an armchair (5,5) CNT under the 
influence of Anderson disorder. Dashed line: g (E) of the clean system displaying the van Hove 
singularities. Solid lines: g {E) from Eq. (11) and £ioc obtained from it via Eqs. (2) and (9). Data 
with errorbars: values obtained numerically by averaging over 100 samples of length 1000 to 5000^uc. 



We model CNTs based on an orthogonal first-nearest-neighbor tight-binding approximation [3] 
with the hopping parameter 70 = 2.66 eV, which accurately reproduces the Fermi velocity of armchair 
CNTs and gives reasonable agreement with the correct band structure for several low energy bands. 

For zigzag GNRs, this approximation is insufficient, as it results in a completely dispersion-free 
edge state while precise calculations show that this state in fact has a small but non-vanishing band 
width [18, 9]. Likewise, this approximation predicts every third armchair GNR to be metallic, while 
precise calculations predict the opening a small gap [19, 20]. The true electronic structure of a GNR 
would have to take into account the relaxed electronic structure at the edges [21]. As it turns out, a 
third-nearest neighbor parameterization [22] of an unrelaxed GNR gives a good approximation of the 
bands near the Fermi energy and will thus be used here. 

Anderson disorder is defined as uncorrelated potential fluctations [10]: 

{EiEj) ~ (Ei) (Ej) =(5jjCre.> (3) 

usually in terms of a uniform random distribution Ei — Eq G [—W/2,W/2], which has a standard 
deviation a^ = W/\/l2. For perturbatively weak disorder, however, the exact shape of the distribution 
is irrelevant. Numerically, we use a Gaussian distribution for practical reasons. Experimentally 
measured values of Id are of the order of 10 /zm in CNT of 1.5 nm diameter [6], which translates to 
an equivalent model disorder of strength a^ ^ 0.05 eV. The values used in Figs. 1 and 2 are chosen 
significantly larger to enhance the visibility of the qualitative effects. 

The diffusive transmission through a quantum wire with a disordered section of length L, 
embedded in an otherwise disorder-free, infinite quantum wire, is defined as 

with the elastic mean free path ^ci- For armchair CNTs at the Fermi energy, this length £c\ was first 
derived by White and Torodov [4] . We have previously given a generalized derivation that is valid for 
arbitrary energies [23] , which can be further generalized to also cover the case of incquivalent atoms 
in the unit cell: 

Within a statistical ensemble P (W) of the perturbation Hamiltonian W, the scattering rate 
between two bands s and d at a fixed energy E is given by the Fermi golden rule 

27r 



giE,d) dWPiW)\{E,d\W\E,s) 



(5) 



which depends on the partial DOS of the destination band d. In the atomic basis \i), the disorder 
Hamiltonian is diagonal and uncorrelated [Eq. (3)], allowing for the direct evaluation of the integral 
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Figure 2. Density of states g{E), number of channels N^h, elastic mean free path £^i and localization 
length £iac of an armchair GNR of width Na = 20 under the influence of Anderson disorder 
(iTj: = 0.2 eV). Dashed lines: g{E) and A^^h of the clean system and lengths obtained from these 
via Eqs. (2) and (9). Solid lines: g (E) from Eq. (11), A^ch from Eq. (13) and lengths obtained from 
these. Data with errorbars: values obtained numerically by averaging over ~ 180 samples of length 
2000 to 20000£uc. 



as 



rE^,^,= '^-^9iE,d)\{E,d\^)fal\{^\E,s) 



(6) 



The elastic mean free path £^^ within a band s is related to the scattering rate r^ ^ from this band via 
its electron velocity which is proportional to the inverse of the DOS as 



''(?-') 



V d J 



(7) 



Based on the definition of £ei via the diffusive transmission in Eq. (4), the elastic mean free path 
^ei of a multichannel quantum wire is found by inverse averaging 



r 



i-E(a- 



A^ch 



(8) 



At this point, the sums over the bands s and d can be reduced, introducing the LDOS on individual 
orbitals gi (E) and resulting in the final expression 

ioi^iucN^iJTT'f^ialgUE))] (9) 



with the length of the unit cell ^u^. and the sum running over all orbitals i within one unit cell. In 
this form, the expression can be applied to arbitrary quantum wires, including GNRs, where it also 
covers the special case of edge disorder by making a^. dependent on the orbital number i. 
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Neglecting multiple scattering, the elastic mean free path 1^1 and the diffusive transmission Tdiff 
are defined in term of the LDOS gi (E) of the disorder-free system. Likewise, iVch is defined by the 
leads, where it follows an exact integer step function. Near band edges, this diffusive transmission 
is discontinuous, as can be confirmed numerically to arbitrary precision, computing it as the sample 
average (T) of the transmission of many disorder configurations [23] . 

The diffusion coefficient, as it can be obtained via the time evolution of a wave packet within 
a long, disordered quantum wire, also allows the extraction of the elastic mean free path [24]. In 
this case, however, it does not depend on the DOS not of the clean, but of the disordered system. 
Near band edges, the DOS depends on the disorder strength non-perturbatively, causing the van- Hove 
singularities to broaden and to shift [25]. This effect has to be taken into account when describing 
diffusion near band edges or around the edge state in GNRs. 

The density of states (DOS) of a quantum wire under the influence of Anderson disorder can 
be obtained via an algorithm based on diagrammatic perturbation theory that takes into account 
localization effects by including multiple scattering [25]. Dropping the crossing diagrams within the 
noncrossing approximation (NCA) [26], allows to write the self energy S (E) to all orders as a recursive 
expression, which can then be iterated numerically until self-consistency is reached. Though the 
applicability of the NCA is not obvious, it be justified by comparing the contribution of various terms 
at low orders [25]. 

The original formulation of this approach is restricted to the special case of CNTs where all 
atoms are equivalent through symmetry and the self energy takes the same value for all atoms. It 
can, however, be generalized to arbitrarily structured quantum wires using matrix notation. The self 
energy E (E) caused by the disorder is a diagonal matrix obeying the recursive relation 

[j:{E)]^^^=S,,al\{E + iO+-no-^{E))-'] . (10) 

For a periodic system, the self energy has the same periodicity as the Hamiltonian. The block- 
tridiagonal matrix (E + iO"*" — TYq — S {E)) can therefore be inverted numerically using a highly 
convergent renormalization-decimation algorithm [27, 23], allowing us to go beyond the energy range 
near the Fermi energy, where the special band structure allows analytic inversion. 

Starting with S = 0, each numerical iteration of this recursive relation is equivalent to one 
additional perturbative order. Typically, convergence is achieved after less then ten iterations, except 
for energies near a van Hove singularity, where hundreds of iterations may be necessary. This clearly 
indicates that low-order perturbation theory breaks down near band edges. 

The LDOS of each orbital i in the unit cell can now be obtained directly from the imaginary part 
of the Green function 

g^E) ^ --lm\{E + iO+- Ho- ^iE)y'] . (11) 

Figs. 1 and 2 show this quantity in direct comparison with the numerically exact value obtained by 
sample averaging. The slight deviation visible at the flanks of the van Hove singularites is caused 
by the NCA [25]. The elastic mean free path £oi based on the DOS of a disordered system, as it is 
displayed in Fig. 2 is no longer a purely perturbative quantity, but it takes into account the scattering 
into localized states present at any given energy. 

The number of channels A^ch in Eqs. (2) and (9) is another quantity that has to be reconsidered 
in the vicinity of band edges. In a disorder-free quantum wire, A^ch is an integer valued step function 
that has discontinuities at band edges. The numerically obtained localization length, on the other 
hand, does not display any discontinuities, so the discontinuities of Net have to be somehow smoothed 
out by the disorder. 

In any periodic system, the number of channels is equivalent to the transmission through a cross 
section. For a Hermitian Hamiltonian, this quantity will always take on integer values. Adding a 
complex self energy may, however, result in non-integer values and smooth out the discontinuities 
of the transmission function. Using the self-energy obtained from Eq. (10) allows thus to define an 
equivalent of the number of channels A^ch for an Anderson disordered quantum wire. 
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Figure 3. Analytically obtained quantities for an armchair edge GNR of width Nz = 41 with 
Anderson disorder of various strengths. Dashed lines: Density of states g {E) of the clean system. 
Solid lines: g (E) and localization length £iqc obtained analytically. Left: homogeneous disorder over 
all atoms across the GNR. Right: disorder concentrated on the outermost atoms along both edges. 
The disorder strength a^^ per atom is scaled such that the square total Tua'^. is the same for both 
types of disorder. 
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The definition of the transmission through a cross section of a quantum wire in the presence of 
a complex self-energy requires some care. Unlike a Hermitian periodic wire, where the transmission 
through a cross section is equal to that through a finite-length section, the imaginary part of the 
self-energy would act as dissipative term in a conducting region and make the result dependent on 
its length. To actually obtain the transmission through a cross section, we should rather split the 
periodic quantum wire in only two parts - similar to the setup of a tunneling junction - writing 

Wl \ / r 
Uk ) yr^ 

with two semi-infinite parts TYl f^nd Tip and the hopping matrix T that connects both parts. Let us 
denote the Green functions of the two isolated semi- infinite leads by Qi^ and Qy{,. The corresponding 
spectral functions are Ai^m, — i^L/R ~ ^^l/r ^^^ ^^"^ Green functions of the complete system 
g-c/a. — (^jt; -^ [Q+ _ 7Y) . If V were a small perturbation like in a tunneling junction its matrix 
elements (^|V|r) would directly give the transition amplitude from an eigenmode TYr jr) = E\r) on 
the righthand side to one TYl |0 = ^ 10 o^i the lefthand side. In a periodic system, however, V is of 
the same magnitude as Tio, so we need to consider the full expansion {l\V + V5''V|r), generalizing the 
tunneling-transmission tr (^lV^IrV) to the nonperturbative expression [23]: 

T = tr [Al (V + Vg'V) Ar (V + Vg^V)] . (13) 

Note that in the case of a finite self-energy S all along the periodic wire, this expression differs 
substantially from more commonly used approaches that depend on the definition of a finite-length 
conductor between two leads. Numerically, we can exploit the finite support of V by computing only 
the contacted regions of Gi, and Qr [27, 23]. 

The localization length .^loc of armchair CNTs at Ep can be derived analytically directly from £e\ 
and iVch [28]. For general energies, £ioc has so far only been accessible by numerical means [29, 30]. 
Based on our expressions for the DOS and the number of channels in a homogeneously disordered 
quantum wire, it is now possible to obtain an approximate value for ^loc in the whole energy range. As 
can be seen in Fig. 1, the analytical value for ^loc agrees with the numerically exact results fairly well. 
The remaining deviation can be understood in view of the oversimplified nature of Eq. (2), which is 
based on a model with a fixed number of equivalent channels with a single value Id for all scattering 
processes. 

Armchair edge GNRs can be understood as unrolled zigzag CNTs and are indeed physically 
very comparable. Similar to semi-metallic (3A'^, 0) nanotubes, which are metallic in zone-folding 
approximation and develop a small gap due to their curvature, ribbons of width Na ~ 3Af — 1 
(counted in rows of carbon dimers) have a tiny gap that is caused by edge effects only. In Fig. 2 
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one can see that the DOS is smoothed out and the van Hove singularities shifted in the same way 
as it was observed in CNTs. Near E-p, this leads to a narrowing of the gap which is well reproduced 
by the data based on Eq. (11). The states near the gap exhibit a short localization length, as it is 
expected from the fact that they are caused by the disorder within the forbidden energy range of the 
clean system. The analytically obtained localization length coincides with the numerical value to a 
procesion comparable to the case of CNTs. Again, the true value is slightly underestimated when 
several inequivalent channels are present at the same energy. 

Reducing the strength of the disorder to realistic values leads to a significant reduction of the 
smoothing effect on van Hove singularieties. Besides this effect, the localization lengths simply scale 
with 1/cr^ wherever band edges have no influence, as can be seen in Fig. 3. Assuming the amount 
of intrinsic disorder in GNRs to be similar to that measured in CNTs at a^ ~ 0.05 eV, we find a 
localization length of (,\oc ^ 13 /^m in armchair GNRs of 5nm width within the plateau of only two 
transmitting channels at sufficient distance to the gap. 

The edges of a GNR must be considered to be especially sensitive to disorder due to the lower 
mechanical rigidity and higher chemical activity. We therefore studied the effect of pure edge disorder 
as depicted in the right panels of Fig. 3. As can be immediately seen, the overall effect of the disorder 
scales with the square total Sicr^. over the disorder strength of all individual atoms. When scaling 
the disorder strength in this way, the effects of edge disorder are identical to the case of homogeneous 
disorder except for the band edges, where the van Hove singularities are smeared out stronger in the 
case of pure edge disorder. 

In zigzag edge GNRs the edge states at the Fermi energy add a considerable complication. 
Numerical studies have been done before on these systems using both Anderson-type disorder of 
short or long range [31, 32] or edge defects [33, 34, 35]. We found that the DOS obtained by Eq. (11) 
shows a considerable deviation from the true value obtained numerically. Possibly, this effect is related 
to the highly anomalous behavior found in two-dimensional sheets of graphene [36] . Trying to obtain 
€ioc via Eq. (2) is bound to fail even more seriously as this approach completely neglects the different 
mean free paths and the spatial separation between the edge channels and the bulk channel. In fact, 
though this approach correctly describes a very short localization length around the edge states, it 
failes to reproduce the correct scaling with the ribbon width. 

To conclude, we have demonstrated that our generalized expression for the density of states in 
Anderson-disordered quantum wires in combination with a highly convergent numerical evaluation 
scheme does reproduce to good precision the true data obtained by performing a numerical sample- 
average. The method allows to efficiently explore of the full energy range and can be applied to 
arbitrary quantum wires, including CNTs and armchair-edge GNRs. We further demonstrated how 
to obtain the localization length of those systems from the density of states in the full energy range, 
including the vicinity of van Hove singularities which are notoriously hard to access by analytical 
means. Again, the resulting data shows good agreement with statistically obtained values. The case 
of zigzag-edge GNRs states has to be excluded due to the extreme inequivalence of the different 
channels close to the Fermi energy and remains a problem for the future. 
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